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A RAPID METHOD FOR OPTIMIZATION 


OF THE ROCKET PROPULSION SYSTEM FOR 
SINGLE -STAGE -TO-ORBIT VEHICLES 

Charles H. Eldred and Sabina V. Gordon 
Langley Research Center 

SUMMARY 

A rapid analytical method for optimization of rocket propulsion systems is presented 
for a vertical take-off, horizontal landing, single -stage -to -orbit launch vehicle. This 
method utilizes trade-offs between system characteristics affecting flight performance 
and engine system mass. The performance results obtained from a point -mass trajec- 
tory optimization program are combined with a linearized sizing program to establish 
vehicle sizing trends caused by propulsion system variations. Although this method is 
applicable and adaptable to a wide range of flight vehicle systems, the example presented 
IS restricted to a single class of vehicle and propulsion systems. The linearized sizing 
technique was developed for this class of systems. The resulting sizing program is 
small enough to be readily adaptable to a desk top programmable calculator and, because 
of its explicit solution, is fast and precise. The specific example treated in this paper is 
the optimization of nozzle expansion ratio and lift-off thrust-to-weight ratio to achieve 
either minimum gross mass or minimum dry mass. Assumed propulsion system char- 
acteristics are high chamber pressure, liquid oxygen and liquid hydrogen propellants, 
conventional bell nozzles, and the same fixed nozzle expansion ratio for all engines on a 
vehicle. Optimization results show a strong trend toward vehicles having low engine 
mass rather than high propulsive preformance, with optimum values of expansion ratio in 
the range of 40 to 45 and optimum lift-off thrust-to-weight ratios in the range of 1.20 to 
1.35. 


INTRODUCTION 

Single -stage-to-orbit (SSTO) systems are being studied as potential launch vehicle 
concepts for the time period following space shuttle. On the basis of projected advance- 
ments in technologies and possible changes in mission requirements, such systems are 
envisioned either as shuttle replacements or in complementary mission roles. If the shut- 
tle system is assumed to have a 15 -year operational life, a new replacement system will 



be required around 1995. With 8 to 9 years allowed for design and development lead time, 
the period available for technology preparation is only slightly more than 10 years. 

Recent studies at the Langley Research Center have been directed toward defining 
the necessary and high-yield areas of technology preparation. An integral part of this 
process is the definition and preliminary design of attractive system concepts, which are 
then used to measure the effects of technology variations. One of the most crucial tech- 
nologies for such systems is that of propulsion, which is the principal subject of this 
paper. 

Some method of total system synthesis is required to compare the various propul- 
sion system concepts and to identify many of the trade-offs within a given propulsion con- 
cept. The extent of synthesis detail utilized in such studies varies over a wide range, the 
most detailed methods requiring large computer resources. The method described in this 
paper emphasizes a simplified synthesis which can reduce the resource requirements for 
such propulsion system trade studies. 

The key to this simplified synthesis is a linearized sizing analysis. Martin, in ref- 
erence 1, uses a linearized sizing analysis to optimize the phasing points of multiphase 
SSTO propulsion systems and to compare various types of propulsion concepts on the 
basis of ideal velocity. The implicit assumption is that the effect of propulsion system 
variations on the ascent trajectory is small in comparison with the variations due to other 
factors. In the present study, a similar type of sizing analysis is used in conjunction with 
actual trajectory performance for analyzing trade-offs between engine systems mass and 
flight performance for a selected engine concept. The sizing analysis is configured for a 
vertical take-off, horizontal landing (VTOHL) SSTO with all rocket propulsion. Either 
gross mass or dry mass may be used as a figure of merit. The fact that the interface 
between the sizing program and the trajectory program is simple allows a different tra- 
jectory program to be substituted if desired. 

The example presented herein is the optimization of engine nozzle expansion ratio 
and lift-off thrust -to -weight ratio to achieve either minimum gross mass or minimum dry 
mass. The propulsion system uses rocket engines operating at high chamber pressures 
with liquid hydrogen and liquid oxygen propellants and conventional bell nozzles, all noz- 
zles on a given vehicle having the same fixed expansion ratio. There are other potentially 
attractive rocket propulsion schemes for VTOHL SSTO application, such as the use of a 
mixture of fixed expansion ratios and dual-position nozzles; however, these are more 
complex than the example scheme, both from considerations of vehicle integration and 
operation and from the standpoint of propulsion system optimization. 
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SYMBOLS 


'sd 


residual propellant mass coefficient, ratio of residual propellant mass to 

*^^p,sd 


ascent propellent mass. 


m 


p,a 


'SV 


reserve propellant mass coefficient, ratio of reserve propellant mass to 


ascent propellant mass, 


m 


p,sv 


m 


P,a 


main engine system mass coefficient, ratio of installed engine system weight 
to thrust, — ^ 


glide -reuse systems mass coefficient, ratio of glide -reuse systems mass to 


mp 

entry mass, — see equation (10) 

rrip 


standard free-fall acceleration, 9.80665 m/s 


engine specific impulse, -7=-, s 


K 


tankage systems mass coefficient, ratio of tankage systems mass to total 


mam propellant mass. 


™p,a *^p,sd + ™p,sv 


m 


engine propellant mass flow rate, kg/s 


m. 


vehicle dry mass, mj + mj^ + mg + mQ, kg 


mg 


engine systems mass, kg, see equation (9) 


m^ 


vehicle entry mass, kg, see equation (13) 


nif 


fixed mass, kg, see equation (7) 


mr 


glide -reuse systems mass, kg, see equation (10) 
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mg vehicle gross mass, kg, see equation (11) 

mj^ tankage systems mass, kg, see equation (8) 

mpay payload mass, kg 

mp 3 ^ mass of ascent propellants, including in-flight losses, kg 

mp o mass of orbital maneuvering propellants, kg 

mp j. mass of reaction control system (RCS) propellants for entry, kg 


’^p,sd mass of main propellant residuals, kg 

”^p,sv mass of main propellant performance reserves, kg 


N 

q 

T 

a 

ra 


vehicle system thrust -to-weight ratio, 


T 

mgg 


free -stream dynamic pressure. Pa 
rocket engine thrust, N 
angle of attack, deg 

ascent propellant mass ratio, ratio of gross mass to mass at ascent trajectory 

nig 

main engine cut-off, 2 

mg - mp^a 


ro 


rr 


orbital maneuvering propellant mass ratio, ratio of initial orbital mass to 
entry mass. 


me + mp^o 


me 


entry RCS propellant mass ratio, ratio of entry mass to landing mass. 


me 


mg - mp r 
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AV 


ideal velocity change, m/s 


e engine nozzle expansion ratio 

Subscripts: 

a main ascent propellant 

o orbit maneuvering system (OMS) propellant 

r entry reaction control system (RCS) propellant 

si sea level or lift-off conditions 

vac vacuum 

Abbreviations: 

LH2 liquid hydrogen 

LINSIZ linear sizing computer program 

LOX liquid oxygen 

ODIN optimal design integration, computerized design integration system 

OMS orbital maneuvering system 

POST program to optimize simulated trajectories, trajectory computer program 

RCS reaction control system 

SSME space shuttle main engine 

SSTO single-stage-to-orbit 

TPS thermal protection system 

VTOHL vertical take-off, horizontal landing 
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METHOD OF ANALYSIS 


Overview 

Propulsion system optimization for a vehicle flight system should include such con- 
siderations as (1) engine performance and mass characteristics; (2) propellant and tankage 
requirements; (3) flight trajectory and performance; (4) integration of the propulsion sys- 
tem with the vehicle (installation, aerodynamic effects, center of gravity); and (5) overall 
system synthesis or sizing. For cases in which the range of propulsion system variables 
under consideration is restricted in some way, a simplified treatment of one or more of 
these factors may be appropriate. The method of analysis presented herein applies sim- 
plifications to (4) and (5) above. Propulsion system variations are assumed not to impact 
configuration features such as vehicle shape, aerodynamic coefficients, vehicle balance 
(center of gravity), etc. Also, system sizing changes are assumed to be adequately rep- 
resented by a combination of linear relationships. However, for both assumptions, this 
does not preclude the need for a reference vehicle design which adequately defines and 
accommodates all the above considerations. 

A principal feature of this method of analysis is the use of a simplified sizing pro- 
gram which incorporates considerations (1), (2), and (5); has the form to permit an explicit 
solution; and is small enough to be run on a desk top programmable calculator. Calibra- 
tion of this sizing program against the characteristics of a more detailed design assures 
proper sizing behavior. 

Ascent trajectory performance, which is assumed to be a driving sizing require- 
ment, is computed with a generalized point-mass trajectory optimization program. Tra- 
jectory performance is determined for various combinations of the propulsion system 
characteristics. Because the effects of integrating the propulsion system with the vehicle 
are assumed not to vary the overall aerodynamic characteristics, it is possible to com- 
pute all trajectories by use of the characteristics of a single vehicle (excluding the pro- 
pulsion system variations of interest). These trajectory results are then used as inputs 
for the sizing program. In this way, the combined trajectory flight performance and 
vehicle sizing effects of propulsion system variations are evaluated. 

Sizing Analysis 

A linearized sizing program was used to determine vehicle sizing for each combina- 
tion of propulsion system characteristics evaluated for trajectory performance. This 
program, designated LINSIZ, was implemented on a desk top programmable calculator in 
about 275 program steps. This linearized method is similar to the approach used in ref- 
erence 1; however, more details in the propulsion area are provided in the present study. 
This sizing approach simplifies the overall vehicle system to 11 mass categories: 
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1 - Mam engine ascent propellant 

2 - In -flight propellant losses 

3 - Ascent propellant reserves 

4 - Ascent propellant residuals 

5 - OMS (and RCS) orbital propellants 

6 - RCS entry propellant 

7 - Payload 

8 - Fixed items 

9 - Tankage systems 

10 - Engine systems 

11 - Glide-reuse structures and systems (wings, TPS, landing gear, etc.) 

Mass groupings that are used in this study are defined as follows: 

Dry mass Categories 8, 9, 10, and 11 

Entry mass Dry mass plus categories 4, 6, and 7 

OMS initial mass Entry mass plus category 5 

Burnout mass Entry mass plus categories 3 and 5 

Gross mass Burnout mass plus categories 1 and 2 (all categories) 

Categories 1 to 7 are self-evident; however, categories 8 to 11 require some expla- 
nation. Vehicle fixed mass includes such items as crew, crew provisions (cabin and 
enviromental control), payload bay structure, payload provisions, and avionics, which are 
independent of propulsion system variations and which are basic to the orbital mission of 
the vehicle. Tankage systems include basic tanks and associated structure and plumbing, 
which are sized by the total ascent propellant requirements. Engine systems include 
engines, plumbing, pumps, gimbals, thrust structure, etc., which are sized by the thrust 
level and propellant flow rate requirements. The glide -reuse category includes wings, 
tail, TPS, landing gear, etc., which are sized by the entry mass. 

For each mass category, a linear mass relationship appropriate to a VTOHL SSTO 
vehicle is used: 
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Ascent propellant, including in-fhght losses, 


- 1 

”^P,a f 

^ a 


Reserves and residuals 

^p,sv + ^p,sd ' (^sv ^sd)*^p,a 
OMS propellant 


mp,o = (^o - 1)^6 


where 


Tn - 6 


'^^o/o^s.o 


RCS entry propellant 


_ Tr - 1 

J^p,r p— me 


( 1 ) 


( 2 ) 


(3) 


(4) 


(5) 


where 


Tj. = 

is a useful approximation. 
Fixed mass 


( 6 ) 


mf = Fj + F2mpay 


(7) 


where Fj is a fixed mass constant (in kilograms), F2 is a fixed mass coefficient 
(fraction of payload mass), and F2mp^y reflects the mass of payload accommodations 
to allow the scaling of fixed mass with payload variations. (For the type of enclosed bay 
with full-length doors assumed in this study, F2 appears to have a value of about 0.6.) 
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However, since payload is not varied in this study, fixed mass is simply m^, 
tion (7) IS not used. 

Tankage systems 

“ ^(^^p,a + ™p,sd *^^p,sv) 

Engine systems 

’^E ~ ^sl^sl*^g 
Glide -reuse systems 

mQ = Gmg 

The definition of gross mass gives the following expression: 

»^^g = + mp^sv + nip,sd + mp^o + nip^j. + mpay + mf + m^ + mg + m^ 

Substitution of equations (2) and (3) gives 

*^g “ ^p,a(^ ^sv) + mp,a^sd + (^o " l)iT^e "*■ *’^p,r + *^pay + r^f + *^E ' 

From the definition of entry mass 

n^e = mp.a^sd + ■"p,r + ™pay + mj + m^ + 
equation (12) can be simplified to 

nig = nip,a(l + ‘^sv) + Torng 

Substitution of equations (5) and (10) into equation (13) gives 

me = T (mp,aCgd + mpay + nif + m^ + m^) 

1 r ~ ri 


and equa- 


( 8 ) 


(9) 


( 10 ) 


( 11 ) 


niQ (12) 


(13) 


(14) 


(15) 
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Then substituting equation (15) into equation (14) gives 


mg = + ^gv) + (mp^aCgd + ^p^y + + nij^ + m^) 


( 16 ) 


1 - -i- G 


Also, substitution of equations (8) and (9) now gives 


nig - nip^a 


1 I o , ^o^sd ^0^ (I I c I C ) 

1 + Cgv + ;; ; + ;; — — (1 + t^gv + <-gc}) 


1 1 - 


- G 


+ 



(mpay + + ^sl^simg) 

- G 


(17) 


And finally, substituting equation (1) and solving for the gross- mass gives 


f TT + m 


f)- 


‘pay 1 / p , 

1 _irjLi 


- G 


mcr = 


1 - 


^o^slNg^ 

- 1 

1 - — i - G 


1 + Cgy + 


^ o*^sd ^o^(l + ^sv + Cgcj) 


1 - 


- G 




(18) 


Equation (18) is the key equation to the explicit sizing solution. Once the gross 
mass is determined, all the other mass characteristics can be established by using the 
following sequence: 

Ascent propellant, including in-flight losses, by using equation (1) 

Reserves and residuals by using equation (2) 

Entry mass by rearranging equation (14) to give 


nif 


iWg - mp^a(^ + ^sv) 


(19) 
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OMS propellant by using equation (3) 


RCS propellant by using equation (5) 

Tankage systems by using equation (8) 

Engine systems by using equation (9) 

Glide -reuse systems by using equation (10) 

Potentially, any of the inputs to equation (18) could be used as design variables for 
a sizing study; however, in applying this sizing analysis to the problem of propulsion sys- 
tem optimization, the inputs fall into three general categories: 


(1) Design assumptions and requirements 


Payload 

OMS requirements 
RCS requirements 
Propellant coefficients 


•"pay 

Pq or AVo and Ig^o 
Pj. or AVj. and Ig^r 
^sv- 


(2) Common mass factors 

Fixed mass ' ~ / 

Glide -reuse coefficient 
Tankage coefficient 

(3) Design variables 

Engine factors 

Thrust-to-weight ratio 
Ascent propellant mass ratio 


mf 

G 

K 


^sl> ^s,a 
Nsl 


After the trajectory performance is determined for an assumed combination of e 
and Ng^, the resulting ascent mass ratio Pg, along with the appropriate engine mass fac- 
tor Egi, is evaluated by the sizing program. The resulting gross mass and dry mass 
are utilized as the overall indicators of propulsion system optimality. 

Within the limitations of the linearized sizing, validity of the results rests primarily 
with the proper selection of values for the three dry mass coefficients E, G, and K 
and the fixed mass mj. In order to establish these coefficient values, a point-design 
vehicle system is needed for which there are detailed mass information and sizing trend 
data. It IS desirable to duplicate both the absolute sizing level and the sizing trends such 
as those associated with variations in mass ratio. 
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Trajectory Analysis 


The purpose of the trajectory analysis is to identify the variations in ascent perfor- 
mance, as measured by the ascent mass ratio which result from variations in char- 
acteristics of the rocket propulsion system. To do this, optimized ascent trajectories 
are computed for each propulsion system configuration by use of a generalized program 
called POST (program to optimize simulated trajectories) described in reference 2. 

POST IS a discrete parameter targeting and optimization program having the capa- 
bility to target and optimize point-mass trajectories while satisfying equality and inequal- 
ity constraints. Because of the generality of the simulation and optimization routines, 
POST can be used to solve a wide variety of problems for a powered or unpowered vehicle 
operating near a rotating oblate planet. These problems can include ascent, reentry, and 
orbital transfer trajectories. In this study, POST is utilized for powered ascent trajec- 
tories which satisfy various in-flight constraints, meet orbital insertion requirements, 
and optimize the vehicle flight controls in order to minimize the ascent propellant require- 
ment. Typical in -flight constraints include limits on dynamic pressure q, product of 
dynamic pressure and angle of attack qof, and vehicle acceleration. Insertion conditions 
may include altitude, flight-path angle, and velocity. 

Interface requirements between the sizing program and the trajectory analysis are 
such that another trajectory program could be substituted for POST if desired. Vehicle 
simulation inputs such as aerodynamics, planform loading, in-flight constraints, and con- 
trol modes must be representative of the vehicle being synthesized in the sizing analysis. 
These inputs are not varied within the trajectory analysis. The inputs that are varied 
between the different trajectory cases are the propulsion system characteristics. The 
sizing -trajectory interface requires that the propulsion system weights used in the sizing 
be consistent with the flight performance characteristics. The principal output of the tra- 
jectory analysis is the ascent propellant mass ratio Fa, which is used as a sizing input. 

Propulsion system characteristics are input as total vacuum thrust, total mass flow 
rate, and total nozzle exit area. Effective thrust due to back pressure losses is calculated 
in the program. In addition to main propellant flows, in-flight propellant losses (due to 
such things as boiloff and auxiliary power generation) can be simulated by the addition of 
a zero-thrust, zero-exit-area "engine." The in-flight losses, which are then included in 
the ascent propellant requirement, have been taken into account in this manner. 

EXAMPLE 

The example presented in this study is the optimization of the rocket engine nozzle 
expansion ratio e and lift-off thrust -to -weight ratio Ng^ for a VTOHL SSTO vehicle 
using identical multiple engines with conventional bell nozzles. All nozzles on a given 
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vehicle have the same fixed expansion ratio. Optimization is on the basis of either mini- 
mum gross mass or minimum dry mass. A Langley Research Center (LaRC) in-house 
preliminary design provided a point-design vehicle for calibration of the linearized sizing 
model and for the trajectory analyses. This point -design vehicle was based on mission 
requirements similar to the space shuttle and on assumed 15 -year advancements beyond 
the space shuttle level in two basic technology areas, structures and propulsion. Trajec- 
tory and sizing analyses were performed for a matrix of e and Ngj combinations, with 
e ranging from 30 to 80 and Ngj ranging from 1.15 to 1.50. 

Design Requirements 

Basic space shuttle design requirements were utilized extensively for this study 
(ref. 3). Although it is unlikely that a post-shuttle vehicle would be designed for the same 
mission needs as the shuttle, the shuttle offers the best defined set of requirements avail- 
able at this time. The sizing mission used was the 29 500-kg payload launched from the 
Kennedy Space Center (KSC) into a 28.5*^ inclined orbit. A 4.6 -m -diameter by 18.3 -m- 
long payload bay is provided. The orbital maneuvering system provided AV = 229 m/s, 
which was sufficient for near -Earth maneuvers and entry deorbit. The OMS provides for 
all maneuvers beyond the 93 -km by 185-km insertion orbit provided by the main propul- 
sion system. The entry cross range requirement is established at the shuttle level of 
2037 km. A reaction control system is provided for attitude control during the initial 
phases of reentry. 


Technology Assumptions 

The assumed development schedule for the SSTO has the advantage of a 15 -year 
advancement in technology beyond the level for the shuttle. Significant gams should result 
from existing programs in the areas of structures and propulsion. 

Structure s. - The information presented in reference 4 indicates that a 25-percent 
structural mass reduction from current technology is possible for a vehicle developed in 
the period from 1985 to 1990. Consequently, a 25 -percent structural mass reduction from 
a space shuttle reference has been assumed in the areas of the body, wing, tail, tanks, and 
landing gear. No mass reductions are assumed for the propulsion systems or for the 
supporting subsystems. 

Propulsion . - The engine characteristics used in this study are based on space shut- 
tle main engine technology, which incorporates a high chamber pressure (20 MPa) and 
LOX-LH2 propellants. A LOX-LH2 mixture ratio of 7:1, as proposed in reference 5, is 
assumed instead of the SSME value of 6:1. This variation appears to offer advantages 
over the SSME for SSTO application because the vehicle sizing advantages due to the 
increased propellant bulk density (about a 9-percent increase) more than offset the effects 
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of the somewhat lower specific impluse (about a 1.5-percent decrease). However, the 
overall differences between these two mixture ratios are not very large. 

Engine characteristics used in this study are shown in figures 1 and 2. The varia- 
tion of specific impulse with expansion ratio e is shown in figure 1. This figure is 
based on the reference point for LOX-LH2 presented in reference 5 (a mixture ratio of 
7:1 and e = 200) and on characteristics of minimum -surface -area nozzles for variations 
in e presented in reference 6. The curve for a mixture ratio of 6:1 with the SSME 
located at e = 77.5 is shown for comparison. Figure 2 shows the variation of the engine 
system mass coefficient E with the expansion ratio for thrust levels at both vacuum and 
sea level. The sea level mass coefficient Egj is the one used in the sizing program. 
These mass characteristics were derived from weight reports and trade studies for the 
space shuttle and reflect the SSME operating at maximum power. 

Point-Design Vehicle 

A point-design vehicle derived from LaRC in-house studies of an SSTO system was 
used for the trajectory evaluation and also to establish coefficients for the linearized siz- 
ing program. This design was generated in Langley Research Center studies which uti- 
lized the ODIN system (ref. 7). Vehicle mass properties and other characteristics are 
summarized in table 1 and the physical layout is shown in figure 3. For this vehicle, e 
was 80 and Ngi was 1.347. The ascent mass ratio F^ determined from trajectory 
analysis was 7.3342. The gross mass was 2.15 Gg, and the reference theoretical wing 
area was 1087 m^. Design requirements and assumptions for the point -design vehicle 
were consistent with those previously discussed, with two minor exceptions. First, the 
OMS mass ratio Fq was insufficient in the point -design vehicle to give the desired AV 

' I 

of 229 m/s, and second, the value of Is assumed for the main engines was slightly higher 
than that assumed for this study, as shown in figure 1. This resulted in a somewhat lower 
ascent mass ratio F^ than would be required for the engines with the lower value of Ig. 
Corrections for both these factors are incorporated in the sizing analysis. 

Sizing Analysis 

Tables 2 and 3 illustrate the use of the point -design vehicle to calibrate the linear 
sizing program. The mass categories of the point-design vehicle were redistributed 
according to the linear sizing categories. The fluid categories were all well defined, how- 
ever, the four categories included in dry mass required some judgment. Assignment of 
the major vehicle components were straightforward, with two possible exceptions: Basic 

structure and Growth. Seventy-five percent of the basic structure (item 3.1) was assigned 
to fixed mass to account for payload accommodations, which include the payload bay, 
structural supports, and doors. Growth (item 18.0) was proportioned among the four dry- 
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mass categories based on the mass level without growth. For many of the smaller sub- 
systems whose categories were not obvious, category assignments were based on corre- 
lations of sizing trends with available trend data. In this current example (and probably 
typical in general), the best sizing trend correlation between LINSIZ and ODIN was 
achieved by placement of these smaller subsystems in the fixed mass category. The sum- 
mation of masses under each category as shown in table 2 provides the information needed 
to establish constants and coefficients for the sizing program. On the basis of this mass 
distribution, the following input values for the calibration vehicle were established: 

Ngi = 1.347 

= 7.3342 (e - 80, Is, a, vac = 453 s) 

Fq = 1.03667 

Fr = 1.0168 

Csv = 0.002 
Cgd= 0.007 
J^pay " 29 500 kg 
m{ = 32 600 kg 


Esi = 0.0248 


K = 0.0240 


G = 0.296 


The resulting calibration vehicle from the LINSIZ program is compared with the 
point -design vehicle in table 3. All categories agree to within 1 percent. Figure 4 com- 
pares the effect on gross mass of varying the ascent mass ratio as predicted by LINSIZ 
and ODIN. The agreement is within 1 percent in the gross mass range of interest for this 
study. However, significant errors are generated farther from the calibration point. The 
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LINSIZ gross mass is higher than the ODIN gross mass at mass ratios above the calibra- 
tion point and is lower below the calibration point. As discussed in the preceding para- 
graph, manipulation of the fixed mass category allows a degree of control on these sizing 
trends, with increases in fixed mass tending to decrease the sensitivity of gross mass to 
changes in mass ratio. 

The baseline vehicle for this sizing study is based on the point-design and calibra- 
tion vehicles but incorporates the corrections for the OMS velocity change AV and spe- 
cific impulse Is of the main engines, as noted earlier. The following changes were 
made to the LINSIZ inputs; 

Ngi - 1.350 


= 7.4536 5 (e = 80, Is,vac = 449.35) 

Tq = 1.05207 (a V o = 229 m/s, Is,o = 460 s) 

The resulting vehicle is about 11 percent heavier than the calibration vehicle and is 
shown in the last column of table 3, The only sizing inputs varied from the baseline vehi- 
cle for the analysis of e and Ng^ effects are Ngj, F^, and Egj. 

Trajectory Analysis 

A typical ascent trajectory is shown in figure 5. The trajectory simulation was set 
up to minimize the ascent propellant requirements while satisfying in-flight constraints 
and orbital insertion requirements. Controls used are pitch angle for the first 30 s and 
the angle of attack a thereafter. These are implemented as linear functions of time in 
eight different segments over the ascent trajectory. The in-f light constraints are maxi- 
mum ascent acceleration, maximum dynamic pressure q, and maximum qcv. Ascent 
acceleration is limited to 3g by continuous engine throttling after the 3g level is reached. 
Dynamic pressure is limited to 48 kPa as an inequality constraint. The no-wind value 
of qo! is limited to 120 kPa-deg by controlling the angle of attack in the region of high 
dynamic pressure. (This qa value was found to be optimal in ref. 8.) Aerodynamic 
coefficients were obtained from a representative Phase B space shuttle orbiter. Orbital 
insertion requirements, which are based on the perigee of a 93- by 185 -km orbit, are 
implemented as equality constraints for altitude and flight -path angle. Velocity is the tra- 
jectory termination parameter. 

In addition to main propellant flows, in-flight propellant losses (due to such things 
as boiloff and auxiliary power generation) were simulated. This amounted to a total loss 
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of 0,17 percent of the gross mass expended at a constant rate over the time of the ascent 
burn. Ascent propellant requirements, which are described by the main propellant mass 
ratio Fa, include these in -flight losses. 

POST iterates the trajectory until convergence criteria which indicate that a satis- 
factory solution has been found are met. Convergence criteria include both constraint 
tolerances and measures of optimality. Typical tolerances on converged solutions are 
2 m on insertion altitude and 0.004° on insertion flight-path angle. 

RESULTS AND DISCUSSION 

The method outlined in this paper for the optimization of SSTO rocket propulsion 
systems is easy and rapid to use once the basic trajectory simulation and sizing program 
constants are set up. The only input changes required between cases for the trajectory 
program are derived from the selection of expansion ratio e and lift-off thrust-to -weight 
ratio Ngp These changes are input as total vacuum thrust, total propellant flow rate, 
and total nozzle exit area. The input changes required for the sizing program are simply 
(1) lift-off thrust -to -weight ratio Ngj, as used in the matching trajectory case; (2) engine 
mass coefficient Egi, corresponding to the selected expansion ratio; and (3) ascent mass 
ratio Fa, as determined by the trajectory. 

Results of the ascent trajectory analysis for the example are shown in figure 6. A 
matrix of cases was evaluated with the expansion ratio ranging from 30 to 80 and the lift- 
off thrust-to-weight ratio varying from 1.15 to 1.50. (The expansion ratio of 80 repre- 
sents an approximate limit for stable nozzle flow at sea level conditions.) The perform- 
ance map of figure 6 utilizes sliding scales for both e and Ng^ in order to have scales 
for both variables. Significant reductions in propellant requirements are possible at the 
larger expansion ratios and at the higher thrust-to-weight ratios. However, increases in 
both expansion ratio and thrust-to-weight ratio are at the expense of increasing engine 
mass. 

Results of the sizing analysis are shown in figures 7 and 8. Figure 7 is a map of 
the dry-mass change relative to the baseline vehicle (e = 80, Ng^ = 1.35). Minimum dry 
mass is achieved at e = 40 and Ng^ = 1.20, a 15.9-percent reduction from the baseline 
vehicle. Details of this vehicle are presented in tables 4 and 5. A comparison of fig- 
ures 6 and 7 shows nearly a complete inversion of results between trajectory performance 
and dry mass. The highest trajectory performance (lowest ascent propellant mass ratio) 
results in the heaviest dry -mass vehicle, whereas the minimum dry mass occurs in the 
region of lowest trajectory performance. It is apparent that the engine mass penalty 
required to achieve the low ascent propellant mass ratios more than offsets the benefits 
of the low mass ratios. 
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Figure 8 is a map of the gross -mass change relative to the baseline vehicle. In this 
case, minimum gross mass is achieved at e = 45 and Ngj = 1.35, a 9.6-percent reduc- 
tion from the baseline vehicle. Details of this vehicle are presented in tables 4 and 5. 

This vehicle provides somewhat higher trajectory performance than the minimum-dry- 
mass vehicle; however, as before, the highest trajectory performance results in the heavi- 
est vehicle. The optimum occurs at a fairly low expansion ratio and a moderate thrust - 
to-weight ratio; therefore, a strong trade-off between trajectory performance and engine 
mass is indicated. 

Comparison of total propulsion system mass in tables 4 and 5 (that is, the sum of 
tankage systems and engine systems)-indicates that the propulsion system mass is on the 
order of 50 percent of the dry mass and that it therefore has both strong direct and resiz- 
ing effects. The propulsion system mass is reduced from 51.7 percent of the dry mass 
for vehicle 1 to 48.4 percent for vehicle 2 and to 49.1 percent for vehicle 3. As a result 
of resizing, the propulsion system mass of vehicle 1 is reduced 21 percent for vehicle 2 
and 17.5 percent for vehicle 3. Another important measure of vehicle optimality is pay- 
load fraction, shown in table 5, based on either dry mass or gross mass. The payload 
fraction based on vehicle dry mass increased from 11.7 percent for vehicle 1 to 14.0 per- 
cent for vehicle 2 and 13.5 percent for vehicle 3. 

Results of this example indicate the strong influence of trade-offs between engine 
system mass and flight trajectory performance for this class of vehicles. These results 
apply only to the VTOHL SSTO vehicle utilizing the assumed propulsion scheme (all 
engines on a vehicle have the same fixed expansion ratio). The optimum vehicles tend to 
accept significant flight performance penalties due to low expansion ratios and low lift-off 
thrust -to -weight ratios in order to achieve the benefits of low engine mass coefficients. 
This illustrates the importance of engine mass coefficient E as a technology area for 
this class of rocket propulsion systems. 

These results also indicate the suitability of this method - linearized sizing com- 
bined with optimized trajectory performance - for the optimization of propulsion sys- 
tems when the range of variables is appropriately limited. As illustrated in the example, 
rocket propulsion systems are characterized by conflicting flight performance and vehi- 
cle sizing considerations and therefore require this type of combined performance-sizing 
analysis for optimization. The linearized sizing analysis provides adequate sizing char- 
acterization while also offering the potential for considerable savings in analysis time 
and/or cost over other sizing analyses. However, because of the linearizing assumptions, 
care must be taken to provide calibration with a suitable vehicle system and to limit the 
extent of sizing perturbations. 

There are a number of other propulsion schemes and operational modes of potential 
interest for SSTO application. Other rocket propulsion schemes include mixed expansion 
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ratios, dual -position nozzles, and dual fuels. The basic assumptions behind this method 
of analysis are quite general; therefore, all these propulsion-operational schemes are 
adaptable to this type of analysis - some require only slight modification of the inputs 
while others may require modifications to the sizing equations. The computation of sen- 
sitivity information is another feature of this method associated with the explicit sizing 
solution. Small perturbations to the sizing imputs can be used to derive sizing sensitivi- 
ties in order to show the impact of technology changes. 

CONCLUDING REMARKS 

A method has been developed for the analytical optimization of the rocket propulsion 
system for a vertical take-off, horizontal landing, single-stage-to-orbit launch system. 
This method combines the performance results obtained from a point -mass trajectory 
optimization program with a linearized sizing program to establish vehicle sizing trends. 
This linearized sizing program, which simplifies the overall system to 11 mass categories 
(4 dry mass categories), has an explicit solution form and is small enough to offer sub- 
stantial computational savings over many other sizing programs. 

Results are presented for a vehicle system utilizing rocket engines with high cham- 
ber pressure, liquid oxygen and liquid hydrogen propellants, conventional bell nozzles, and 
the same fixed expansion ratio for all engines on a vehicle. Both expansion ratio and lift- 
off thrust-to-weight ratio were optimized to achieve vehicles with either minimum gross 
mass or minimum dry mass. Results show that propulsion system mass is a dominant 
vehicle sizing factor. The trades between engine system mass and trajectory perform- 
ance show a strong trend toward low engine mass. Optimum values of expansion ratio 
fell in the range of 40 to 45; optimum values of lift-off thrust-to-weight ratio fell in the 
range of 1.20 to 1,35. 

With proper care in its application, both in establishing the required sizing inputs 
and in limiting the range of perturbations, this method should give quite satisfactory 
results. Also, the method is adaptable to many other flight vehicle/propulsion system 
concepts. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
May 19, 1976 
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TABLE 1.- POINT-DESIGN VEHICLE - MASS PROPERTIES 


VTOHL SSTO with LOX-LH2 propellant; reference wing area, 
1087 m2; Ngi = 1.347; e = 80; and Is, vac = 453 s 


Item 

Component category 

Mass, kg 

1.0 

Wing group 

23 717 

2.0 

Tail group 

3 454 

3.0 

Body group (total mass, 62 515 kg) 


.1 

Basic structure 

19 575 

.2 

Thrust structure 

5 538 

.3 

LOX tanks 

12 296 

.4 

Fuel tanks 

25 106 

4.0 

Induced environmental protection 

3 5 882 

5.0 

Landing, docking, recovery 

8 735 

6.0 

Propulsion - ascent 

61 291 

8.0 

Propulsion - auxiliary 

3 825 

9.0 

Prime power 

1 775 

10.0 

Electrical, conversion and distribution 

3 403 

11.0 

Hydraulic, conversion and distribution 

3 171 

12.0 

Surface controls 

3 305 

13.0 

Avionics 

2 021 

14.0 

Environmental control 

1 857 

15.0 

Personnel provisions 

790 

18.0 

Growth 

15 475 


Dry mass 

231 216 

20.0 

Personnel 

705 

21.0 

Cargo 

29 483 

22.0 

RCS reserves 

68 

23.0 

RCS propellant 

4 536 

24.0 

Residuals 

13 064 


Entry mass 

279 072 

25.0 

OMS propellant 

10 224 

26.0 

Reserve fluids 

3 725 

27.0 

Mam propellant (total mass, 1 855 998) 


.1 

In-flight losses 

3 725 

.2 

Ascent 

1 852 273 


Gross mass 

2 149 019 
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TABLE 2.- POINT -DESIGN VEHICLE - BREAKDOWN TO LINEAR SIZING CATEGORIES 
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TABLE 3.- MASS PROPERTIES OF POINT -DESIGN VEHICLE AND RESULTING LINSIZ DESIGNED VEHICLES 
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TABLE 4.- SUMMARY OF MASS PROPERTIES FOR BASELINE VEHICLE AND TWO OPTIMIZED CONFIGURATIONS 
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TABLE 5.- MASS DISTRIBUTION FOR BASELINE VEHICLE AND TWO OPTIMIZED CONFIGURATIONS 
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Relative velocity, km/ s 


Dynamic pressure, kPa 


Ascent propellant mass ratio. 
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Figure 6.- Effects of expansion ratio and lift-off thrust-to-weight ratio 
on ascent trajectory performance. 


29 



30 40 50 60 70 80 

Expansion ratio, e 


Lilt 

1.50 1.35 1.25 1.15 

Sea level thrusHo-weight ratio, 


Sliding scales 


Figure 7.- Effects of expansion ratio and lift-off thrust -to-weight ratio 

on SSTO vehicle dry mass. 
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